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ABSTRACT: In this paper, we introduce an approach to increase density of field-effect transistors 
framework current mode instrumentation amplifier flipped voltage follower mirrors. The approach we 
consider for manufacturing of the inverter in heterostructure with specific configuration. Several required 
areas of the heterostructure should be doped by diffusion or ion implantation. After that dopant and 
radiation defects should by annealed framework optimized scheme. We also consider an approach to 
decrease the value of mismatch-induced stress in the considered heterostructure. We also introduce an 
analytical approach to analyze mass and heat transport in heterostructures during manufacturing of 


integrated circuits with account mismatch-induced stress. 


KEYWORDS: current mode instrumentation amplifier; optimization of manufacturing; analytical 


approach for modelling. 


https://doi.org/10.29294 /IJASE.8.1.2021.2046-2068 ©2021 Mahendrapublications.com, All rights reserved 


INTRODUCTION 

In the present time several actual problems of 
the solid state electronics (such as increasing of 
performance, reliability and density of elements of 
integrated circuits: diodes, field-effect and bipolar 
transistors) are intensively solving [1-6]. To 
increase the performance of these devices it is 
attracted an interest determination of materials 
with higher values of charge carrier’s mobility [7- 
10]. One way to decrease the dimensions of 
elements of integrated circuits to manufacturing 
them in thin film heterostructures [3-5,11]. In this 
case, it is possible to use inhomogeneity of 
heterostructure and necessary optimization of 
doping of electronic materials [12] and devel- 
opment of epitaxial technology to improve these 
materials (including analysis of mismatch induced 
stress) [13-15]. An alternative approaches to 
increase the dimensions of integrated circuits are 
using of laser and microwave types of annealing 
[16-18]. 


In this paper, we introduce an approach to 
manufacture field-effect transistors. The approach 
gives a possibility to decrease their dimensions 
with increasing their density framework a current 
mode instrumentation amplifier flipped voltage 
follower mirrors. We also consider the possibility 
to decrease mismatch-induced stress to decrease 
quantity of defects, generated due to the stress. 


We consider a heterostructure, which consist of a 
substrate and an epitaxial layer (see Fig. 1). We 
also consider a buffer layer between the substrate 
and the epitaxial layer. The epitaxial layer includes 
into itself several sections, which were 
manufactured by using another material. These 
sections have been doped by diffusion or ion 
implantation to manufacture the required types of 
conductivity (p or n). These areas became sources, 
drains and gates (see Fig. 1). After doping it is 
required annealing of dopant and/or radiation 
defects. 


The main aim of the present article is to analysis 
of redistribution of dopant and radiation defects 
to determine the conditions, which correspond to 
decreasing of elements of the considered filter and 
at the same time to increase their density. At the 
same time, we also consider a possibility to 
decrease mismatch-induced stress. 


Method of solution 


To solve our aim we determine and analyzed 
spatio-temporal distribution of concentration of 
dopant in the considered heterostructure. We 
determine the distribution by solving the second 
Fick's law in the following form [1,19-23] 
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x) Oa =a 
Fig. 1a. Structure of the considered amplifier [24] 
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Fig. 1b. Heterostructure with a substrate, epitaxial layers and buffer layer (view from side) 
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x=L, 
Here C(x,y,z,t) is the spatio-temporal distribution 
of concentration of dopant; Q is the atomic volume 
of dopant; Vs is the symbol of surficial gradient; 


i 
(C(x, y,z,t)dz is the surficial concentration 
0 


of dopant on interface between layers of 
heterostructure (in this situation we assume, that 
Z-axis is perpendicular to interface between layers 
of heterostructure); su(x,y,z,t) is the chemical 
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Here Dz (x,y,z,T) and Dis (x,y,z,T) are the spatial 
(due to accounting all layers of heterostructure) 
and temperature (due to Arrhenius law) 
dependences of dopant diffusion coefficients; T is 
the temperature of annealing; P (x,,z,T) is the 
limit of solubility of dopant; parameter y depends 
on properties of materials and could be integer in 
the following interval ye[1,3] [21]; V(x,y,z,t) is the 
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potential due to the presence of mismatch-induced 
stress; D and Ds are the coefficients of volumetric 
and surficial diffusions. Values of dopant 
diffusions coefficients depends on properties of 
materials of heterostructure, speed of heating and 
cooling of materials during annealing and spatio- 
temporal distribution of concentration of dopant. 
Dependences of dopant diffusions coefficients on 
parameters could be approximated by the 
following relations [21-23] 
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spatio-temporal distribution of concentration of 
radiation vacancies; V is the equilibrium 
distribution of vacancies. The concentrational 
dependence of dopant diffusion coefficient has 
been described in details [21]. Spatio-temporal 
distributions of concentration of point radiation 
defects have been determined by solving the 
following system of equations [19,22,23] 
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Here I (xy,z,t) is the spatio-temporal 
distribution of concentration of radiation 
interstitials; [* is the equilibrium distribution of 
interstitials; Di(x,y,z,T), Dv(x,y,z,T), Dis(xy, z,T), 
Dys(x,y,Z,T) are the coefficients of volumetric and 
surficial diffusions of interstitials and vacancies, 
respectively; terms V2(x,,z,t) and /[(x,y,z,t) 
correspond to generation of divacancies and 
diinterstitials, respectively (see, for example, [23] 


and appropriate references in this book); 
kiv(xy,Z,T), ku(xy,z,T) and kvv(x,y,z,T) are the 
parameters of recombination of point radiation 
defects and generation of their complexes. 


Spatio-temporal distributions of divacancies ®y 
(x,y,z,t) and diinterstitials ® (x,y,z, t) could be 
determined by solving the following system of 
equations [19,22,23] 
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Here Da(x,y,z,T), Dov(x,y,z,T), Dass (x,y,Z,T) and Devs(x,y,z,T) are the coefficients of volumetric and surficial 
diffusions of complexes of radiation defects; ki(x,y,z,T) and kv(x,y,z,T) are the parameters of decay of 
complexes of radiation defects. 


Chemical potential sm in Eq.(1) could be determine by the following relation [19] 
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tensor; ui, uj are the components u:(x,y,z,t), u(xy,z,t) and u:(x,y,z,t) of the displacement vector 
Uu (, y,Z; t); Xi, Xj are the coordinate x, y, z. The Eq. (3) could be transform to the following form 


where E(z) is the Young modulus, oy is the stress tensor; uj, = is the deformation 
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where oa is Poisson coefficient; & = (as-azz)/az is 
the mismatch parameter; ds, az, are lattice 
distances of the substrate and the epitaxial layer; 
K is the modulus of uniform compression; / is the 
coefficient of thermal expansion; T; is the 


equilibrium temperature, which coincide (for our 
case) with room temperature. Components of 
displacement vector could be obtained by solution 
of the following equations [25] 
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We determine spatio-temporal distributions of corrections [26]. Previously we transform the 
concentrations of dopant and radiation defects by Eqs.(1), (3) and (5) to the following form with 
solving the Eqs.(1), (3) and (5) framework account initial distributions of the considered 
standard method of averaging of function concentrations 
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and radiation defects in right sides of Eqs. (1a), for the first-order approximations of the required 
(3a) and (5a), on their not yet known average concentrations in the following form 
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Integration of the left and right sides of the Eqs. (1b), (3b) and (5b) on time gives us possibility to obtain 
relations for above approximation in the final form 


om V (x, ¥.2, Vek 
C(x, Vrzt)= M% QA | Dslr Hah |b a Dee oe ] x 


ane ae Ot ee 
Vsti (tsy.267)) 14+ 9 Nd ch aye | Dyy (x, 95257) b+ 9 
x sosare)} “PGs Gen) sil 927) "Prxy.2T)| 


2 
x OV u(x, y, ea ieee (ny ce aun a oe | dt+ f.(x,y¥.Z) Gd 


04D 0o+tD 
AY et)= My 2O5 [TEV seo Wet )d T+ O20-— JE V tales vet) de 
oy (a y> z)— a fk, (x, y, 2,0 )d 0 @,Qy [py (x, y> z,T) dt (3c) 
0 0 
01D O24) 
Vibe ys zot)= Oy 22 |e Vstilen Weve hd t+ Oy ZO [FEV seule woz) det 
+ fy (19.2) ai [ky yO Ys2T)d tO, Oy Sk, y(% ZT) dr 


6 1Dos 6 tDoys 
D(x, y,z,t)=a,, ZA — PVM, y,Z,7) dt+Q_2QM— “V.UW\x y,zZ,7) dx 
a y ) 10, a! kT 5M ( y ) ae! kT an » ) 


XQhy, Z+ fy, (x, yZ)+ +{k, (x, y,z,T)I (x, y,z,7)d r+{k,, (x, y,z,T )I’(x, y,2z,7)dt (50) 
0 0 


0 1 Do,s 0 + Dos 
DO, (x, y,z,t)=A, ZA— “Vv x, y,z,7)dtT+Q— “Vv x,y,Z,T)dTX 
iv y ) loy rae kT stl y ) a kT fal y ) 


Xie, Z+ Ie (x, y, z)+ [k, (x, y, zT)V (x, y, z,T)d T ‘ie. (x, y, z,T V(x, y, z,T)d 4 
0 0 


Pankratov 


International Journal of Advanced Science and Engineering www.mahendrapublications.com 


Int. J. Adv. Sci. Eng. Vol.8 No.1 2046-2068 (2021) 2055 E-ISSN: 2349 5359; P-ISSN: 2454-9967 


We determine average values of the first-order approximations of concentrations of dopant and radiation 
defects by the following standard relation [26] 
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We determine approximations of the second and concentrations of dopant and radiation defects we 
higher orders of concentrations of dopant and replace the required concentrations in the Eqs. 
radiation defects framework standard iterative (1c), (3c), (5c) on the following sum anptp ni 
procedure of method of averaging of function (x,y,z,t). The replacement leads to the following 
corrections [26]. Framework this procedure to transformation of the appropriate equations 
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Integration of the left and the right sides of Eqs. (1d), (3d) and (5d) gives us possibility to obtain relations 
for the required concentrations in the final form 
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Average values of the second-order approximations of required approximations by using the following 
standard relation [26] 
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Further, we determine the solutions of Eqs.(8), i.e. function corrections we replace the required 
components of displacement vector. To determine functions in the right sides of the equations by 
the first-order approximations of the considered their not yet known average values a. The 
components framework method of averaging of substitution leads to the following result 
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Integration of the left and the right sides of the above relations on time t leads to the following result 
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Approximations of the second and higher orders of components of displacement vector could be 


determined by using standard replacement of the required components on the following sums aj+ui(x,y,z,t) 
[26]. The replacement leads to the following result 
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Integration of the left and right sides of the above relations on time t leads to the following result 
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To support, in this paper we also determine 
concentration of dopant, concentrations of 


radiation defects and components of displacement 
vector by using the second-order approximation 
framework method of averaging of function 
corrections. This approximation is usually enough 
good approximation to make qualitative analysis 
and to obtain some quantitative results. All 
obtained results have been checked by 
comparison with results of numerical simulations. 


DISCUSSION 

In this section, we analyzed the dynamics of 
redistributions of dopant and radiation defects 
during annealing and under influence of 
mismatch-induced stress. Typical distributions of 
concentrations of dopant in heterostructures are 
presented on Figs. 2 and 3 for diffusion and ion 
types of doping, respectively. These distributions 
have been calculated for the case, when value of 
dopant diffusion coefficient in the epitaxial layer is 
larger, than in the substrate. The figures show, 
that inhomogeneity of heterostructure gives us 
possibility to increase sharpness of p-n- junctions. 
At the same time one can find increasing 
homogeneity of dopant distribution in doped part 
of epitaxial layer. Increasing of sharpness of p-n- 
junction gives us possibility to decrease switching 
time. The second effect leads to decreasing local 
heating of materials during functioning of p-n- 
junction or decreasing of dimensions of the p-n- 
junction for fixed maximal value of local overheat. 
However framework this approach _ of 
manufacturing of bipolar transistor it is necessary 
to optimize annealing of dopant and/or radiation 
defects. Reason of this optimization is following. If 
annealing time is small, the dopant did not achieve 


any __ interfaces between materials of 
heterostructure. In this situation one cannot find 
any modifications of distribution of concentration 
of dopant. If annealing time is large, distribution of 
concentration of dopant is too homogenous. We 
optimize annealing time framework recently 
introduces approach [15,25-32]. Framework this 
criterion we approximate real distribution of 
concentration of dopant by step-wise function (see 
Figs. 4 and 5). Farther we determine optimal 
values of annealing time by minimization of the 
following mean-squared error 
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Fig.3. Distributions of concentration of implanted 
dopant in heterostructure from Fig. 1 in direction, 
which is perpendicular to interface between 


epitaxial layer substrate. Curves 1 and 3 
corresponds to annealing time © = 
0.0048(Lx2+Ly2+L22)/Do. Curves 2 and 4 
corresponds to annealing time © = 
0.0057(Lx2+Ly2+Lz2)/Do. Curves 1 and 2 


correspond to homogenous sample. Curves 3 and 
4 corresponds to heterostructure under condition, 
when value of dopant diffusion coefficient in 
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epitaxial layer is larger, than value of dopant diffusion coefficient in substrate 


C(x,0) 


Fig.4. Spatial distributions of dopant in heterostructure after dopant infusion. Curve 1 is idealized 
distribution of dopant. Curves 2-4 are real distributions of dopant for different values of annealing time. 
Increasing of number of curve corresponds to increasing of annealing time 
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Fig. 5. Spatial distributions of dopant in heterostructure after ion implantation. Curve 1 is idealized 
distribution of dopant. Curves 2-4 are real distributions of dopant for different values of annealing time. 
Increasing of number of curve corresponds to increasing of annealing time 
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Fig.6. Dependences of dimensionless optimal annealing time for doping by diffusion, which have been 
obtained by minimization of mean-squared error, on several parameters. Curve 1 is the dependence of 
dimensionless optimal annealing time on the relation a/L and €= y= 0 for equal to each other values of 
dopant diffusion coefficient in all parts of heterostructure. Curve 2 is the dependence of dimensionless 
optimal annealing time on value of parameter ¢ for a/L=1/2 and €= y= 0. Curve 3 is the dependence of 
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dimensionless optimal annealing time on value of parameter € for a/L=1/2 and ¢ = y= 0. Curve 4 is the 
dependence of dimensionless optimal annealing time on value of parameter yfor a/L=1/2 and e= €=0 
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Fig.7. Dependences of dimensionless optimal annealing time for doping by ion implantation, which have 
been obtained by minimization of mean-squared error, on several parameters. Curve 1 is the dependence of 
dimensionless optimal annealing time on the relation a/L and €= y= 0 for equal to each other values of 
dopant diffusion coefficient in all parts of heterostructure. Curve 2 is the dependence of dimensionless 
optimal annealing time on value of parameter ¢ for a/L=1/2 and €= y= 0. Curve 3 is the dependence of 
dimensionless optimal annealing time on value of parameter € for a/L=1/2 and ¢ = y= 0. Curve 4 is the 
dependence of dimensionless optimal annealing time on value of parameter yfor a/L=1/2 and e= €=0 
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Fig. 8. Normalized dependences of component uz of displacement vector on coordinate z for nonporous 
(curve 1) and porous (curve 2) epitaxial layers 
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Fig. 9. Normalized distributions of charge carrier mobility in the considered heterostructure. Curve 1 
corresponds to the heterostructure, which has been considered in Fig. 1. Curve 2 correspond to a 
homogenous material with averaged parameters of heterostructure from Fig. 1 
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where yw (x,y,z) is the approximation function. 
Dependences of optimal values of annealing time 
on parameters are presented on Figs. 6 and 7 for 
diffusion and ion types of doping, respectively. It 
should be note that it is necessary to anneal 
radiation defects after ion implantation. One could 
find spreading of concentration of distribution of 
dopant during this annealing. In the ideal case, 
distribution of dopant achieves appropriate 
interfaces between materials of heterostructure 
during annealing of radiation defects. If dopant did 
not achieve any interfaces during annealing of 
radiation defects, it is practicably to additionally 
anneal the dopant. In this situation optimal value 
of additional annealing time of implanted dopant 
is smaller, than annealing time of infused dopant. 


Further, we have analyzed influence of 
relaxation of mechanical stress on distribution of 
dopant in doped areas of heterostructure. Under 
following condition a< 0 one can_ find 
compression of distribution of concentration of 
dopant near interface between materials of 
heterostructure. Contrary (at s>0) one can find 
spreading of distribution of concentration of 
dopant in this area. This changing of distribution 
of concentration of dopant could be at least 
partially compensated by using laser annealing 
[29]. This type of annealing gives us possibility to 
accelerate diffusion of dopant and other processes 
in annealed area due to inhomogenous 
distribution of temperature and Arrhenius law. 
Accounting relaxation of mismatch-induced stress 
in heterostructure could leads to changing of 
optimal values of annealing time. Mismatch- 
induced stress could be used to increase density of 
elements of integrated circuits. On the other hand 
could leads to generation dislocations of the 
discrepancy. Fig. 8 shows distributions of 
component of displacement vector, which is 
perpendicular to interface between layers of 
heterostructure. 


CONCLUSION 


In this paper, we modeled redistribution of 
infused and implanted dopants with account 
relaxation mismatch-induced_ stress during 
manufacturing field-effect hetero transistors 
framework a current mode _ instrumentation 
amplifier flipped voltage follower mirrors. We 
formulate recommendations for optimization of 
annealing to decrease dimensions of transistors 


2067 E-ISSN: 2349 5359; P-ISSN: 2454-9967 


and to increase their density. We formulate 
recommendations to decrease mismatch-induced 
stress. Analytical approach to model diffusion and 
ion types of doping with account concurrent 
changing of parameters in space and time has 
been introduced. At the same time the approach 
gives us possibility to take into account 
nonlinearity of considered processes. 
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